UG-06-02 



Dynamics of Generalized Assisted Inflation 



Jelle Hartong, Andre Ploegh, Thomas Van Riet and Dennis B. Westra 

Centre for Theoretical Physics, University of Groningen, Nijenborgh 4, 
9747 AG Groningen, The Netherlands 

■ j. hartong, a. r. ploegh, t. van. riet, d.b.westra@rug.nl 

o ■ 

(N 

^ • Abstract 
^. 

, We study the dynamics of muhiple scalar fields and a barotropic fluid in an 

fS| ' FLRW-universe. The scalar potential is a sum of exponentials. All critical points are 
constructed and these include scaling and de Sitter solutions. A stability analysis of 

1^ [ the critical points is performed for generalized assisted inflation, which is an extension 

I of assisted inflation where the flelds mutually interact. Effects in generalized assisted 

■ inflation which differ from assisted inflation are emphasized. One such a difference 
, is that an (inflationary) attractor can exist if some of the exponential terms in the 

O ' potential are negative. 
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1 Introduction 



Scalar fields can violate the strong energy condition, which is a necessary condition for 
inflation and present-day acceleration. An alternative way to violate the strong energy 
condition is pure vacuum energy, that is, a positive cosmological constant A. Observations 
do not exclude either possibility, but scalar fields provide more freedom to evade typical 
fine-tuning problems posed by a cosmological constant such as the smallness of A or the 
cosmic coincidence problem. On the other hand it is argued that scalar fields generate new 
problems not posed by a cosmological constant [1]. 

In this paper we investigate scalar field models and we confine ourselves to models where 
the scalar fields have a canonical kinetic term and the potential is a sum of exponential 
terms where the exponent is a linear combination of the scalars. The reason for this choice is 
twofold. First of all, these models allow to find exact solutions, which correspond to critical 
points in an autonomous system ^. The second reason is that exponential potentials arise 
in models motivated by string theory such as supergravities obtained from dimensional 
reduction (see for instance [3-7] and references therein) , descriptions of brane interactions 
[8-10], nonperturbative effects and the effective description of string gas cosmology (see 
for instance [11]). 

The simplest example is the single-field model with one positive exponential potential, 
which was shown to have a stable powerlaw inflationary solution when the potential is 
sufficiently flat [12]. Later on this model was extended by including a barotropic fluid, 
negative exponentials, spatial curvature and combinations thereof [13-17]. In all cases it is 
possible to rewrite the equations of motion as an autonomous system. Some critical points 
correspond to so-called scaling solutions. Scaling solutions are defined as solutions where 
different constituents of the universe are in coexistence such that the ratios of their energy 
densities stays constant during evolution. 

Since supergravities typically contain many scalars we are interested in the extension 
of these models to multiple scalars. In reference [18] the authors consider a potential that 
is a sum of exponentials with each exponential containing a different scalar. The main 
result is that if each separate exponential is too steep to drive inflation the system can 
still have an inflationary attr actor provided the number of fields is sufficiently large. This 
effect is known as assisted inflation. In reference [19] this system is extended to include a 
barotropic fluid and spatial curvature. 

In assisted inflation there is no mutual interaction (cross-coupling) between the scalars. 
To understand the effects of mutually interacting scalars reference [20] investigates a model 
with one exponential term containing multiple scalars. But as we will argue this theory 
can be rewritten as a theory without cross-coupling. Mutual interaction of the scalar fields 
can only be present when there are more exponentials terms in the potential. 

As shown in [21] multiple exponential potentials are naturally divided into two classes: 

1. The minimal number of scalar fields which can occur in the potential equals the 
number of exponentials. 

^For a review on dynamical systems in cosmology see [2]. 



2 



2. The minimal number of scalar fields which can occur in the potential is less than the 
number of exponentials. 

A more precise definition can be found in section|21 The first class can have scaling solutions 
and the second class can have scaling and de Sitter solutions. 

Here we shall mostly be interested in class 1. This class is termed generalized assisted 
inflation in [22]. Generalized assisted inflation includes assisted inflation as a special case 
where there is no cross-coupling. In reference [22] some exact solutions of generalized 
assisted inflation were given, which were later shown to correspond to a subset of all 
possible critical points [21]. 

In contrast to assisted inflation, generalized assisted inflation is poorly investigated. 
Therefore we extend the results of [21,22] in two ways. First, by constructing all possible 
critical points of the model containing the most general multiple exponential potential with 
a barotropic fluid and spatial curvature. The second and more important extension is a 
stability analysis of the critical points in generalized assisted inflation. 

The paper is organized as follows. In section 2 the equations of motion are given in a 
flat FLRW-background and are then rewritten as an autonomous system. In section 3 we 
construct all critical points for the most general multiple exponential potential. Section 4 
contains a full stability analysis of all critical points in generalized assisted inflation. In 
section 5 the model is generalized by the inclusion of nonzero spatial curvature and finally 
in section 6 we summarize and conclude. In the appendices we prove some results. 



2 The model 



The scalars are denoted by 
by the following action 



and are assembled in an A^-vector 6. The model is defined 



S 



— 7^- 1(9^0, 



M 



(1) 



where = 8nG with G Newton's constant and {d^cp , d^cp) is shorthand for Xlili d^(pid^(pi 
The potential V{(f)) is a sum of M exponential terms 



M 



(2) 



a=l 



where (a^ , 0) = J2iLi Oiai4>i- There are M vectors with components Uai- The indices 
j, . . . run from 1 to A^ and denote the components of the vectors and a^. The indices 
a,b, . . . run from 1 to M and label the different vectors aa and the constants Aq. Sm is 
the action for a barotropic fluid with density p and pressure P and as equation of state 
P = (7 — l)p. It is assumed that the barotropic fluid respects the strong energy condition, 
which means that 2/3 < 7 < 2. 
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In this paper we make use of linear field redefinitions. If the scalars transform linearly 
as — > (f)' = S(j), where S is an element of GL(]R, A^) then the vectors transform in the 
dual representation aa — > a'a = S~'^aa- This can be seen from the definition of a'^ 



an 



(«1 , 0') • (3) 



Field redefinitions that shift the scalar fields leave the aa invariant, but change the A^. 

From the action we can deduce some properties of this system by looking at transfor- 
mations in scalar space. The kinetic term is invariant under constant shifts and 0{N)- 
rotations of the scalars. These transformations map the multiple exponential potential to 
another multiple exponential potential but with different and aa- Such redefinitions 
do not alter the physics, they just rewrite the equations. Therefore qualitative features 
only depend on 0(A^)-invariant combinations of the Oa-vectors (for example , ab))- By 
shifting the scalars we can always rescale R of the A^ to be ±1, where R is the number of 
independent a-vectors. 

The Ansatz for the metric is that of a flat FLRW-universe and accordingly the scalars 
can only depend on the cosmic time r 

ds'^ = -dr"^ + a{TYdx'^ , 0^ = 0i(r) . (4) 

In section |S1 the analysis is extended to spatially curved FLRW-universes. The equations 
of motion are 



i(0,0) + V(0) + pJ, (5) 

<Pi + ?,H<Pi + diV{(l)) = , (6) 

p + 37i^p = 0, (7) 

where H = d/a is the Hubble constant. The dot denotes differentiation with respect to r. 
The equations (jMZ|) are referred to as the Friedmann equation, the Klein-Gordon equation 
and the continuity equation, respectively. 

One can rewrite equations (jSHZj) as an autonomous system. We define the following 
dimensionless variables 

^. = ^, >; = ^Aaexp[-«:(a.,0)], = ^ • (8) 

If we write = ^ . and Y = J2a the equations of motion become 

n + X^ + Y -1 = 0, (9) 

X; = 3X,(-l + X2 + 2n) + ^1 (10) 

a 

Y; = Ya(^-VQ {aa , X) + ^ g^^j ^ ^^^^ 

= 3^(^2x2 + 7(^]- 1)) , (12) 
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where the prime denotes differentiation with respect to ln(a) ^. 

It can be shown that if f2 + X"^ + Y — 1 = initially then equations p0lll2j) guarantee 
that it is satisfied at all times. Hence, given the correct initial conditions the dynamics is 
described by equations p()lll2|l . 

The number of linearly independent vectors aa is denoted by i?. If i? < one can 
rotate the scalars such that (pR+i, . . . ,4>n no longer appear in the potential {aai = for 
i > R). These scalars are then said to be decoupled or free. 

In the next section we will construct all critical point solutions of the autonomous 
system fll()lll2|l . For that purpose it is useful to consider the matrix A 

Aab = {aa , ab) ■ (13) 

The models are divided into two classes. The first class is defined by i? = M and the 
second class by i? < M. Algebraically the two differ in the following way 

1. R = M detA>0, (14) 

2. R<M detA = 0. (15) 

The first possibility, where A is invertible, is called generalized assisted inflation. 

Assisted inflation is the subclass where A is diagonal. This implies that in assisted 
inflation the aa are perpendicular to each other and that one can choose an orthonormal 
basis in which = aa^ai- In that basis the potential becomes 

M 

V{(j)) = ^Aaexp[-Kaa(j)a\ ■ (16) 

a=l 

It is this particular form of the potential that is referred to as assisted inflation in the 
literature. We want to emphasize that the latter definition is basis-dependent. Potentials 
different from ()lfi|l but with a diagonal matrix A can be brought to the form (jlfjj) through 
an 0(A^)-rotation of the scalar fields. An 0(A^)-invariant definition of assisted inflation is 
that A is a diagonal matrix. 

In any system with multiple fields but a single exponential such as studied in [20], the 
matrix A is trivially diagonal. One can perform a rotation on the scalars such that only 
one scalar appears in the potential and all the others are decoupled. In order to have a 
system whose scalars are mutually interacting one needs at least two exponential terms 
both containing more than one scalar. 

^The use of ln(a) as a time coordinate fails when a = 0. This is no problem for studying critical points 
and their stability because in the neighborhood around a critical point there always exists a region where 
the coordinate is well defined. In reference [16] an explicit example is given where d becomes zero at some 
point. 
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3 The critical points 



Critical points are defined as solutions of the autonomous system for which fl' — X[ — 
— 0. From the acceleration equation, 



(0,0)+7P , (17) 



it follows that in a critical point H/H^ is constant. If the constant differs from zero we 
put H / H"^ = —1/p and then the scale factor becomes 

a(T) = ao(-)^ (18) 
To 

In terms of the dimensionless variables p can be expressed as 

_ 1 

^ ~ 3(^2 + in) ■ ^^^^ 

When H / H"^ = the scale factor is 

a(T)=aoe^% (20) 

and space-time is de Sitter. 

The requirement that — can be satisfied in two ways as can be seen from 



Y; = Ya {aa , X) + QX^ + 37!^) . (21) 



Either = or the second factor on the right hand side equals zero. If we put = by 
hand and then solve for the Xi and the remaining Ya, the critical point is called a nonproper 
critical point. If we put the second factor to zero by hand and then solve for X^ and Ya, the 
critical point is called a proper critical point. The name nonproper is given since a critical 
point with some Ya = has cxo-valued scalar fields. Therefore these critical points are no 
proper solutions to the equations of motion, they are asymptotic descriptions of solutions. 
The proper critical points generically have nonzero Ya and therefore are proper solutions 
to the equations of motion. But in some cases one finds that Fq = for proper critical 
points, although one did not put those Ya to zero by hand. 

Regardless of whether critical points are proper sohitions to the equations of motion, 
they are all equally important in providing information about the orbits. That is, they are 
either repellers, attractors or saddle points. 

3.1 Generalized assisted inflation 
Proper critical points 

First we construct the proper critical points by proceeding as explained in [21] and find 
Ya^2^^Y.^A-%, X, = -JlY,aa^[A-X,, O = 1 - - J^^^U . (22) 

"^^6 P ab P ab 
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The value of p is found by combining equations (011211221) • One gets 



p'-p{2 Y,[A-\, + + Y.'^A-X, = . (23) 



ab ab 

The two possible solutions are 



2J2[A''U and P<A+P = |^- (24) 



ab 



The first solution has vanishing barotropic fluid {Q = 0) and is called the scalar- dominated 
solution. The second solution is a matter- scaling solution [15]; the ratio of the fluid density 
Q to the scalar field density fl^ = k'^p^/SH'^ (with ^ , (f)) -\-V) is finite and constant. 
For the matter-scaling solution the scalar field mimics the fluid resulting in p = ^ , which 
is the same powerlaw for a universe filled with only the barotropic fluid. 
In terms of the scalar fields the solutions read 

(pi = ln|r|+Q. (25) 



Nonproper critical points 

To construct nonproper critical points we put a certain subset of Ya equal to zero by hand 
and then solve for the Xi and the remaining Ya. We denote the elements of the subset of 
Ya that were put to zero with barred indices a, so Fa = 0. 

The solutions are still given by equations (j^^ and (f^^ but now the indices a, b, c run 
over the subset of nonzero Ya and the matrix A has to be replaced by the matrix A{d,b, . . .), 
which is the submatrix of A formed by deleting the d,b, . . . columns and rows of A. 

Similar to the proper critical points, the nonproper critical points are divided into 
matter-scaling solutions and scalar-dominated solutions. If all Ya = the Friedmann 
equation Q gives that + = 1. From ()12|) we notice that either = or = 0. For 
the first case where f2 = 0, the solution is called a kinetic- dominated solution and p = 1/3. 
For the second case where f2 = 1, the solution corresponds to a universe only filled with a 
barotropic fluid and is termed the fluid- dominated solution {pp = 2/87). 

The maximal number of critical points can be seen to be 

Max H critical points = 2 x 2^^ (26) 

The reason is that there are 2*^ different ways of putting Ya equal to zero. For each 
combination there are two possibilities: when not all Ya equal zero there is the scalar- 
dominated and the matter-scaling solution and when all = there is the fluid-dominated 
solution and the kinetic-dominated solution. 

For the scalar-dominated solutions the nonproper critical points correspond to universes 
with an expansion that is always slower than that of the proper critical point. This is due 
to the following lemma. 
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Consider an invertible matrix A with as entries innerpro ducts of vectors Aab = («a , 
Define A{a, b, . . .) as before, then we have the following inequality 

J2i^-'U>Y.l^{a,b,--VU. (27) 

ab cd 

For diagonal A this property is obvious. The proof for general A can be found in appendix 

El 

Existency conditions 

The critical points constructed above do not always exist. For instance, it is clear from 
the definition of the yQ-variables that they must have the same sign as the A^. So, if the 
solution for does not respect that, it has to be discarded^. Of all critical points two 
always exist, the kinetic-dominated and fluid-dominated solution. For the other critical 
points it depends on whether they are matter-scaling or scalar-dominated solutions. 

For the existence of a matter-scaling solution there is an extra condition coming from 
r2 > 0. For the matter-scaling solution the conditions Q > and sgn{Ya) = sgn(Aa) 
become 

E[^"'U<i' JL'^ab^O for A„^0. (28) 

ab ^ 

This shows for instance that for a single negative exponential potential there is no matter- 
scaling solution [16]. 

If the critical point is a scalar-dominated solution then existence is guaranteed if 

{—-3)J2[A'%b^0 for A„^0. (29) 

P<t> 

For the nonproper critical points the same equations hold by changing A to A{a, b . . .). 

To conclude, in generalized assisted inflation we distinguish four kinds of critical points: 
the (non)proper scalar-dominated solutions, the (non)proper matter-scaling solutions, the 
fluid-dominated solution and the kinetic-dominated solution. The kinetic and fluid-dominated 
solutions always exist and for the matter-scaling and the scalar-dominated solutions the 
existency conditions are given in and (jSHI), respectively. 

3.2 The R< M case 

This section is brief since details can be found in [21]. 

An important difference with the previous case where R = M is that these potentials 
can have stationary points diV{(f)) = which correspond to de Sitter solutions (a ~ e"^) 

•^In [4] it is shown that critical points that violate the existency conditions can still play a rol in 
understanding the late time behavior of general solutions. 
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when V > 0. Since R < M the «„- vectors are hnearly dependent and we can always 
choose an independent set of R vectors with a — 1 . . .R. The remaining vectors 
{b — R. . . M) can be expressed as hnear combinations of that independent set 



R 



a=l 

Concerning the proper critical points one has to distinguish between the case of potentials 
with affine coupling'^ and those without. By affinc coupling we mean the situation where 
we can find a set of independent eta's such that the coefficients c^a obey 

X C6„ = 1 for all h . (31) 

a 

These couplings are for instance found in group manifold reductions of pure gravity [4] . If 
the couplings are affine then the proper critical points correspond to powerlaw solutions. 

Since the inverse of the matrix A does not exist when R < M the formulas for the 
solutions differ. We introduce the symmetric matrix B 

M 

Bij = X aaiaaj . (32) 

a=l 

If we assume that a basis rotation is performed whereby the decoupled scalars are removed 
from the potentials we can restrict the index i to run from 1 to R. The decoupled scalars 
obey = . . . = Xn = 0. With this restriction on the index i the matrix B becomes an 
R X R invertible matrix. The powerlaw solutions when the coupling is affine is given by 

^^ = ^E^^^'««^- (33) 
For p there are again two possibilities 

i ja 

Therefore there exist proper scaling solutions (matter-scaling and scalar-dominated) when 
the number of exponential terms exceeds the number of scalars, but only if the coupling is 
affine. When R < M and the couphng is not affine there are no proper scaling solutions, 
but there can still exist nonproper scaling solutions (which are referred to as approximate 
scaling solutions in reference [5]). 

When the couplings are not affine the only possible proper critical point is a de Sitter 
universe. 



^In [21] we called this coupling a convex coupling. Later we found that the name afRne is more 
appropriate. 
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For the construction of the nonproper critical points a subset of the Ya must equal zero. 
Depending on the subset that is put to zero, the resulting equations are that of the R = M 
or R < M-case. So in general nonproper critical points in the R < M case can be de Sitter 
and powerlaw. 

The existence of the critical points is more complicated than for generalized assisted 
inflation since the Fa-variables are no longer independent variables but are nonlinear de- 
pendent. 

4 Stability analysis 

To determine the stability of a critical point one linearizes the equations of motion around 
the critical point. The matrix that appears in the linearized autonomous equations is 
called the stability matrix. Stability of the critical point requires that the real parts of the 
eigenvalues of the stability matrix are negative. 

We will not discuss the stability of critical points for models with R < M but only for 
R = M. The reason is the difficulty that arises from the dependence of the Ya. It is a 
hard problem to write down the independent perturbations around a critical point. There 
exist exceptions where the perturbations are independent when R = M — 1, see [23] for an 
example where R = 1, M = 2. 

Assume (Xj, Ya, Q) is a critical point and consider a perturbation {Xi+6i, Ya+6a, Q+6n)- 
If we substitute the perturbations in equations (fTn|l and (fTT|) and keep the linear terms we 
find 

S^ = il- 3)5. + 3(2 - J)J2X^X,6, - J](|7^i5a " ^faa^^a) , (35) 

S'a = 5^(6(2 - 7)X,- - VQaaj)s, + 6a(^-VQ{aa , X) + ^ - 37^) . (36) 
j 

The Sn perturbations are obtained by varying the Friedmann equation 0. Under this 
assumption it is not an independent perturbation. The equations above can be written 
as a matrix equation 6' = M6. From (j^l36|) one can read off the values of the stability 
matrix M 

Mi,- = [J-3]% + 3(2-7)X,X,, 

M,, = -f X, + ^ati , (37) 
Maj = 6(2 - ^)YaXj - VGaajYa , 

Mab = lab (- V^(«a , X) + ^) - 371^ • 

We now discuss the decoupled scalars 0_r+i, . . . ,4>n- Since Oaii+i = . . . = aaN = it 
follows from equation (|Tn|l that either Xr+i = . . . = = or some Xi^n ^ 0. The 
latter only happens for the kinetic-dominated solution (X^ = 1, = 0). The effect of free 
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scalars at the kinetic-dominated solution is discussed in [4]. For all critical points except 
the kinetic-dominated solution M has only diagonal components in the directions i > R: 

^ M^j ^ , Mij = - 3)Sij for i>R. (38) 

The decoupled scalars decouple in the stability-analysis and do not introduce instabilities 
as long as p > 1/3. 

Wc first investigate the eigenvalues of the stability matrix at the nonproper critical 
points and then at the proper critical points. 

4.1 Stability of nonproper critical points 

The sub matrix Mij is diagonalizable with an -transformation that rotates the scalars. 
The eigenvalues are 

E,^[l-3] , = [i - 3] + 3(2 - 7)X^ (39) 

The multiplicity of i^i is i? — 1 and the multiplicity of E2 is 1. This is sufficient to 
understand the stabihty of the critical points for which all Ya — 0; the kinetic-dominated 
and the fiuid-dominated solution. For these two solutions we then diagonalize the Mij- 
submatrix via an 0(-R)-rotation and we end up with a matrix M which is upper triangular. 
The eigenvalues can then be read off from the diagonal. For the fluid-dominated solution 
the eigenvalues are 3(7/2 — 1) < and 37 > 0. Therefore the fluid-dominated solution is 
always unstable. For the kinetic-dominated solution the eigenvalues are 3(2 — 7) > and 
—^/6{aa , X) + 2/p, the flrst eigenvalue is always positive and hence the kinetic-dominated 
solution is always unstable. 

Assume a certain nonproper critical point has one Ya equal to zero and all other Ya 
nonzero. After renumbering we can choose it to be Yi. The Yi-perturbation decouples 
from all the other perturbations. The reason is that ii Yi — the a — 1 row of M only 
contains its diagonal element Ma=i a=i and therefore the determinant of M — A becomes 

det(M - A) = iMa=ia=i - A)det(M - A) , (40) 

where M denotes the matrix obtained after deleting the a = 1 row and a = 1 column 
from M. The matrix M is the same as the stability matrix of the proper critical point 
belonging to the new system that one obtains by removing the Aiexp[— fi;(ai , 0)]-term 
from the potential. 

This means that the stability of the nonproper critical point with Yi — and Ya^i ^ 
boils down to checking the sign of Ma=\ a=i and to study the stability of the proper critical 
point of the truncated system M. The stability of proper critical points is left for the next 
section and we now focus on the sign of Ma=ia=i- 

Stability requires 

Ma=ia=i = -V6{ai ,X) + 2/p<0. (41) 
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In appendix^ we prove that this expression is equivalent to 

J2[A-'U<0. (42) 

b 

For nonproper critical points with more Y^'s equal to zero, Yi = Y2 = . . . = Yq = we 
find C conditions 

R-C+l R-C+l 

5^ [A(2,...,C)]-,i<0 J2 [A{l,...,C-l)]-,'<0. (43) 

b=l b=l 

Sb[^~^]ab > for all a = 1 ... M then equation cannot be satisfied (see appendix 
IX]) and there does not exist a stable truncation of the Ya. In the case of a diagonal A matrix 
this property is obvious, hence we find that in assisted infiation nonproper critical points 
are always unstable. 

4.2 Stability of proper critical points 

The characteristic polynomial 

If we consider perturbations around a proper critical point a simplification occurs in the 
stability matrix. The Mab submatrix reduces to Mat = —37^0, since now 

("-^) = 4- ^44) 

The eigenvalues of M are found by solving the characteristic polynomial det(M — Al) = 
0. The computation is simplified using an LU-decomposition of M — Al. The characteristic 
polynomial becomes 

det(Mi,- - Xlij) det{Zab) = , (45) 
where we defined the matrix Z 

Zab = Mab - Mab " M;[(M - Xl)-\' . (46) 

The LU-decomposition is not valid when A is an eigenvalue of Mij (-£^1,2) • Solutions to 
det{Z) = are eigenvalues of M if they differ from -Ei,2- With this LU-decomposition one 
has to check separately whether Ei^2 are eigenvalues of M. 

The matrix Zab contains only contracted i-indices and is therefore 0(-R)-invariant. We 
can choose a specific orthonormal basis to compute it. The easiest choice is that basis in 
which Mij is diagonal, where we have 

M;[(M - Xl)-']/M,b = -^^M:M,b + - -^^WaRMu, . (47) 

This expression contains and UaR in the specific basis. One can determine them 
explicitly since by definition of the new basis X = (0, . . . , 0, \fX^). Using this and the fact 
that expression is 0(-R)-invariant we find that aaR = 2/pV6X^ and hence 

Zab = cYa + ^ ^ ^ AabYa - X6ab , (48) 

-C/1 — A 
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with c the following expression^ 



c = -87 + 



E2-X 



9(2 - 7)7^^ + 



37-12 
P 



+ 



p2X2j p'^X^iEi-X) 



(49) 



If one adds all rows to the first row and then subtracts the first column from all the other 
columns the new matrix Z' one obtains has the same determinant as Z and is given by 



Z' 



cY + 



1 



p2 Ei-X 



-A 



Ei-X 



Ei-X 





(^afe>l — Aal) — ^abX 



(50) 



This shows that the determinant factorizes in a c-dependent and a c-independent polyno- 
mial 

det{Z')= (^cY + ^^^-^-\yet(^^^{Aab-Aai)-5at\) witha,6>l. (51) 

The factorization of the characteristic polynomial in a c-dependent and a c-independent 
part allows us to divide the perturbations into fluid perturbations and scalar field pertur- 
bations. The fluid perturbations are those that give rise to the eigenvalues that come from 
the c-dependent factor whereas the scalar field perturbations belong to those eigenvalues 
that come from the c-independent factor. 



Fluid perturbations 

If we put the first factor of ()5H) to zero we get the following eigenvalues for the scalar- 
dominated solution 

Ai = Ei, A2 = -- 37, (52) 
P 

and for the matter-scaling solution 



3 



4 



A± = A = ^(2-7) -l±Wl-87 ^ -'^a,i- ja. _ (53) 
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These eigenvalues have a simple interpretation. Stability of the scalar-dominated solu- 
tion requires 

A2 < ^ p> 2/37 > p^+p. (54) 

For stability it does not matter whether Ei is an eigenvalues or not, because its implication 
for stability {p > 1/3) is less strong then the implications coming from A2 {p > 2/37) when 
we assume that 2/3 < 7 < 2. 

Stability of the matter-scaling solution requires 

ReA±<0 p>2'^[A-%h P4>+p>P4>- (55) 

ab 



^For the special case of i? = 1 the formulae can be taken over if one replaces Ei by E2 everywhere. 
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The existency condition for the matter-scahng solution ()28|) imphes stabihty ()55p which in 
turn imphes instabihty of the scalar-dominated solution. In [24] it is proven that for la- 
grangians that allow for scaling solutions but without cross-coupling, the scalar-dominated 
solution is unstable when the matter-scaling solution is stable. 



Scalar field perturbations 

The remaining eigenvalues follow from the smaller polynomial 

det { ^^^_^ i^ab -Aai)- 6abX) = with a,b>l. (56) 

Determining the solutions of for arbitrary A is a hard problem and we will use another 
technique to determine stability. This technique does not use the (linearized) autonomous 
equations and was developed in the context of assisted inflation [18] and elaborated upon 
in [25]. Below we will extend this to the case of generalized assisted inflation. 

The technique is based upon making a basis rotation in scalar space such that one 
direction corresponds to the critical point solution and the other directions are constructed 
perpendicular to it. Then the idea is to check whether the critical point direction is a 
minimum of the potential with respect to the other direction, if so, the critical point is an 
attractor. 

In order to understand which basis rotation to make notice that the proper critical 
points are given by 

(pi{r) = ^^'P In a , (57) 

K 

irrespective of whether it is a matter-scaling or a scalar-dominated solution. Combinations 
of the form 

^ - ^ , (58) 
Xi Xj 

are constant in time at the critical point. Therefore we define the following new scalar 
combinations 

0: = iV,[^ - |±i] for t<R, (59) 



R R\. v v V V v V V V V v 

A2A3...AK A3A4...A_rAi AiA2...A_r_i 



with Ni and A^^^; normalization constants given by 

-1/2 



NR = X-'lltiXi. (61) 



The directions (p'i^n are all perpendicular to the direction 0'^, but they are not perpen- 
dicular among themselves, hence this is not an orthogonal field redefinition. To establish 
orthogonality we perform a Gram-Schmidt procedure on the directions. The Gram- 
Schmidt procedure takes linear combinations of the to obtain new directions 0^^^ 
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that are mutually perpendicular. It follows from the construction that the are also 
constant at the critical point. We define 0'^ = 0'^. 

The transition of the (pi to the (f)'l is done orthogonally and therefore the system still 
has a canonical kinetic term and a multiple exponential potential 

S = / d'xV^[—n - , dy") - 5^ A, exp[-«: {a'^ , 0")] J • (62) 

a=l 

We focus on the appearance of 0'^ in the potential, so we have to find what a"^ is. 
The scalars are redefined by an orthogonal transformation, 0" = 00. The a-vectors 
transform as dual vectors a'^ = [0~^]'^ Oa which simplifies to a'^ = O Ua for orthogonal 
transformations. In components this means that 

= O/aaj • (63) 

To find a'^ji it then suffices to construct for all i. 

Although we do not know the matrix O since we did not specify the Gram-Schmidt 
procedure, we can find O}^ from (p?j) since 0j 



R-^R 



Thus in the potential the scalar 0'^ appears in the a-th exponent as 

2< 



c^'Lr^r = X-\X , a,) 0'A = . (65) 



The product a"^0'^ is independent of a, so all exponentials have the same coupling to 0'^ 
and we can rewrite the potential as 



M 



V{<P") = e v^.- W{<P';, . . . , 0'^_J , W{<P';, . . . , 0'^„i) =J2^a e~-<<'^"^ , (66) 



a=l 



and (a" , 0") does not contain 0'^; 



R-l 



i=l 

The fact that we were able to rewrite the potential like ()66p extends the results of [25] 
to the case of generalized assisted inflation. 

The function W{(j)i, . . . , 4''ji_i) in ()66|1 contains R—1 scalar fields and has R exponential 
terms. These functions can have a stationary point {dV = 0) and this stationary point 
is unique. Moreover this stationary point corresponds to the critical point, the solution 
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for which the fields 0'/, . . . , take constant values. This is proven by considering the 
Klein-Gordon equations in the new basis 

^"i + 3H(t)'\ + di{e^^ W) = 0. (68) 

Since for the critical point 0", . . . , 0'^_i are constant we see that 

^,MW{cP';,...,<P'^^^)) = 0. (69) 

The critical point is an attractor if and only if this stationary point is a minimum of W. 
Stability requires that the matrix didjW{4>'(, . . . , 0'^_i) (for i,j = 1 . . . ,R — 1) is positive 
definite. If all > then didjW is positive definite since didjW = a^rja with r] a diagonal 
matrix with positive entries, rjaa = Aaexp[— ^(a^ , 0")] > 0. Hence if the proper critical 
point of a potential with all > exists, it is stable. When some A^ < we have proven 
in appendix IbI that the matrix didjW = a^rja is not positive definite if p > 1/3. Hence 
potentials with some negative A^ cannot have stable proper critical points with p > 1/3. 

4.3 Summary of the stability analysis 

We discuss qualitatively the results first for assisted infiation and then for generalized as- 
sisted inflation. We explain how to find the stable critical points (attractors) by eliminating 
all unstable critical points. 

Assisted infiation 

At the end of section H?T] we found that in assisted inflation all nonproper critical points are 
unstable. So if an attractor is present it is either the proper matter-scaling or the proper 
scalar-dominated solution. 

Because of the existency conditions for critical points with diagonal A (f^l^ . the 
proper critical point only exist when all Aq > 0®. This implies that for assisted inflation 
models with some negative A^ there is no attractor, all critical points are unstable. 

In subsection 14.21 we analyzed the eigenvalues coming from the c-dependent factor of 
the characteristic polynomial det(Z) = 0. The stability conditions are p^ > p^+p for the 
scalar-dominated solution and p^ < p^+p for the matter-scaling solution. This means that 
the stable proper critical point is the one that represents the universe with the fastest 
expansion of the two. The nonproper critical points are all unstable and have smaller 
p-values, which implies that the attractor is the critical point with the fastest expansion 
of all the critical points. Furthermore, when the matter-scaling solution exists, it is stable 
and the scalar-dominated solution is unstable. 

In sum, assisted inflation models have an attractor when all exponential terms are 
positive and furthermore the attractor is unique and corresponds to the critical point with 

^We assume that p > 1/3, since scalar-dominated solutions with > 1/3 are unstable because one 
can proof that _Ei = 1/p — 3 is always an eigenvalue of the stability matrix for assisted inflation. 
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the fastest expansion. Thus in assisted inflation the universe dynamically favors a fast 
expansion! 

A fast expansion is favored because of the 'assisting behavior' of the scalar fields. Each 
scalar field gives a positive contribution to the Hubble friction and thus more scalar fields 
lead to more friction and more friction leads to slow rolling fields. Therefore one says that 
the scalars assist each other in helping the universe to expand. The assisting behavior is 
proven for all multi-field lagrangians that allow for scaling solutions and have no cross- 
coupling of the scalars [24]. 

Generalized assisted inflation 

We separately discuss the cases where all terms in the potential are positive and where 
some are negative. 

1. V Aa > 0. If a proper critical point exists, the scalar field perturbations are stable 
since all A^ > 0. As in assisted infiation the fiuid perturbations imply that the 
stable proper critical point is the one that represents the universe with the fastest 
powerlaw, either the matter-scaling (if p^+p > P(f,) or the scalar-dominated solution 
(if > Pp+4,). Since the proper critical point exists we have ^b[^~^]afe > for all 
6^, which implies that truncations are unstable (see discussion below (jlSjl ). Hence 
the nonproper critical points are unstable and the attractor is unique. The cases 
analyzed in [7] are of this type. 

If there exists no proper critical point, attractors have to be nonproper. To find 
stable nonproper critical points, the stability of the truncation ()42p43j) is the only 
condition that has to be fulfilled. Therefore nonproper critical points can be stable 
when the proper critical point does not exist ^. It is not guaranteed that there is 
a unique stable truncation nor if a stable truncation is possible at all. Hence the 
number of attractors is case dependent and can be zero, one or higher. 

2. 3 Afl < . Whereas critical points with negative do not exist in assisted inflation, 
they can exist in generalized assisted inflation but are unstable as we now explain. 
If . . . , are negative at a critical point p than so are A„^, . . . , A„^ in order 
for p to exist. If p is proper we proved at the end of section 4.2 that it must be 
unstable. If p is nonproper i.e. = . . . = Yg^ = then we showed in section 4.1 
that the , . . . , y^p-perturbations decouple from all other perturbations. The latter 
perturbations are the same perturbations as that of a proper critical point p of a 
lagrangian where A^^ , . . . , A^^ are put to zero. This truncated lagrangian contains 
the negative , • • • , A„^ hence p is unstable. Since the perturbations around p 
contain the perturbations around p, p must also be unstable. This proves that 
proper and nonproper critical points with negative Ya are unstable. 

^We use that > 1/3 for scalar-dominated solutions of positive potentials. 

^But the kinetic- and fluid-dominated solution can never be stable as we showed in section 4.1. 
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Assume Ai, . . . , A„ < and all other A^ > 0. If some of the Yi, . . . ,Yn differ from 
zero they are negative in order for the critical point to exist. Therefore we have to 
truncate at least all Yi, . . . ,Yn to find a stable critical point. If such truncations are 
unstable, there is no attractor. Again, it is not guaranteed that there is a unique 
stable truncation nor if a stable truncation is possible at all. 

5 The effects of spatial curvature 

We show that for critical points with k = curvature perturbations do not affect the 
eigenvalues found in section 4. Instead curvature perturbations lead to an additional 
eigenvalue. 

The metric Ansatz is 

ds^ = -dr' + air f (-^ + r'dQ') , (70) 
1 — fcr^ 

where dQ"^ is the metric on the unit 2-sphere and A; = 0, ±1 is the normalized curvature. 
The equations of motion are 



i(0,0) + V(0) + pl -4, (71) 

0, + 3H(j), + diV{(f)) = , (72) 
p + 3-fHp = . (73) 

Using the same dimensionless variables as in section El the equations of motion become 



n + X^ + Y -1 = -^, (74) 



X[ = X,(g - 2) + ^1 ^ aaiYa , (75) 

a 

YI = Y, (- {a, , X) + 2(g + 1)) , (76) 

n' = n(2Q - {3^ - 2))y (77) 



where in the last line a new quantity Q, the deceleration parameter, is introduced. It is 
defined by 

Q^-^ = 2X'-Y + ni^-l). (78) 



For powerlaw solutions a ^ the deceleration parameter is simply Q = 

One might wonder if the system is still autonomous due to the -^-term. Using equations 
()75II77|) one shows that 

(x^ + F + ^]-l)' = 2g(x^ + F + ^]-l). (79) 
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The solution to this equation is 



\X' + Y + Q-l\ = ^, (80) 



with C > 0. If equation (fTIjl is satisfied at the instant r = tq, it is satisfied at all times. 
We conclude that the evolution never changes the value of k and that the dynamics is still 
described by the autonomous system ()75ll77p . 

The critical points can be divided into critical points with = and critical points 
with 7^ 0. The critical points with = are the same as the ones discussed in section IHl 
The critical points with A; 7^ necessarily have p = 1. 

Curvature perturbations 

We define the function /(X, F, Vt) by 

/(X,y,fi) =X2 + F + ^]-l. (81) 

The function / provides us with information about the spatial curvature, for example 
sgn/ = k. Consider perturbations Xj + 5j, Ya + 5a and VL + 5^- They induce a perturbation 
5f which we will call a curvature perturbation. It satisfies the following differential equation 

{5fy = 2Q5f + 2f5Q. (82) 

In section0](/ = 0) we considered perturbations that satisfy 5f = 0. In order to study the 
stability of the critical point against curvature perturbations we need to treat Si, Sa and 
5q as independent perturbations. The stability matrix M can be calculated in exactly the 
same manner as in section 14.21 with the only difference that the resulting matrix has one 
extra row and one extra column. The entries are 



Mnn = 2g + (37 - 2) (fi - 1) , Mnj = SnX, 



] 1 



Man = (37 - 2)Ya , = (§7 - l)^i , 

Maj = -VQYaaaj + SYaXj , Mnb = -2^ , ^ (83) 



Mij = {Q-2) hj + 4X,Xj , = -Xi + ^ lau , 

Ma, = -2Ya + (2(g + 1) - aa^X,) 1„6 . 

Performing an LU-decomposition of the matrix M — Al one shows that 

det(M - Al) = (Mnn - A)det(Mi, - Ali,)det(a,)det(Z,fe) = 0, (84) 
where Cij and Zab are given by 

a. = % - 5:(M., - Al..)~^^^ , (85) 

Za. = Ma, - Xla, - - C-^[M,, - A1.)-(M.„ - (M. - . 
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This decomposition is valid if 



A^Mnn, Xy^E,, Xy^E2, det(a,) ^ 



^6) 



where Ei = Q — 2, E2 = Q — 2 + are the eigenvalues of Mij. Making use of the basis 
in which X = (0, . . . , 0, Vx^) as explained in section |3] we find 



Zab 



4fi 



37 



cYa + 



E2-X Mnn 
3 



A 



Ei-X 



AabYa — Xlab ' 



(87) 
(88) 



where the constant c is 
37-2 



2 + 2n- 



2 {Q + lf 



X2 



A 



5(Q + l)fi 



El- X 
37 



X2 

2 



nn 



X 



+ 8nx 



detCjj E2 - 
2 37-2 



HQ 



nn 



X 



-4X2 
37-2 



nn 



X 



^9) 



The matrix Z is of the same form as in section 14.21 but with a different constant c. In 
section it was shown that for proper critical points the determinant of Zab factorizes in 
a c-dependent and a c-independent polynomial. This result also applies to proper critical 
points with /c 7^ 0. Therefore the only eigenvalues that change are those that are related 
to c. These eigenvalues are determined by 



X = cY + 



6p-2 



El - A 



(90) 



This proves that the effects of curvature only manifest themselves in a subset of the 
eigenvalues of the stability matrix. In this sense the effects of curvature decouple from the 
scalar dynamics. 

Next we discuss the effect of curvature perturbations for the scalar-dominated and 
scaling solutions with A; = as well as the stability of the new critical points with k ^ 0. 
In both cases we only consider proper critical points. For the case of nonproper critical 
points the arguments of section 14.11 apply here as well. 



Stability of critical points with k = 

For the scalar-dominated solution we find the eigenvalue 

2 

Ai = Mnn = 87, 

P 

directly from the expression for M (jHS))- Equation (jnn|) is solved by 

2 



X2 — El , A3 



- 2. 



P 



(91) 



(92) 
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The value of A2 violates the conditions 
stronger constraint than £"1 < 0. 

For the matter-scaling solution equation (j9n|l gives 



however A3 < requires p > 1, which is a 



5(2 

4^ 



7) 



-1 ± 



-1 



ab 



7 



A3 = 37 - 2 , A4 = M, 



nn ■ 



(93) 



Eigenvalue A4 violates the conditions ()86j) . However, the condition A4 < is equivalent to 
A3<0. 

In both cases the eigenvalues Ai^2 are also found in section and the eigenvalue A3 is 
additional. The condition A3 < implies Q < -v^ a > 0. That Q < has to be interpreted 
as stability against curvature perturbations is clear from equation (jHH)- It follows that the 
inflationary critical points are stable against curvature perturbations. The matter-scaling 
solution is never inflationary since the fluid respects the strong energy condition. 



Stability of critical points with /c 7^ 

For the proper critical points with 7^ we have Q = Q = from which it follows Y = 
2X^ = I J2ab ^ab- '^^^ curvature can be found by computing /(X, Y,^l) = 2 J^ab ^ab " ^■ 
The solutions to equation are 



Ai,2 = -1± /8 5^A„-,^-3. (94) 

y ab 

Since f2 = 0, A3 = 2 — 37 is an eigenvalue. Stability requires that 

ab 

Hence the only stable fixed points are the ones that have k = —1. These are the curvature- 
scaling solutions of reference [17]. By a curvature-scaling solution is meant that the energy 
density that effectively describes the effects of spatial curvature is a constant fraction of 
the total energy density. 

If one puts all Ya equal to zero, then all Xi are zero and the critical point is never an 
attractor and represents a Milne patch of Minkowski space-time. 



6 Conclusion 

In this paper we study the stability of critical points in models with multiple exponential 
potentials and multiple fields. The critical points are separated in proper and nonproper 
critical points. The proper critical points are exact solutions to the equations of motion 
and the nonproper solutions are asymptotic descriptions of solutions. A stability analysis 
is carried out for models that belong to the class of generalized assisted inflation. In that 
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class a critical point with A; = is either a kinetic-dominated solution, a fluid-dominated 
solution, a matter-scaling solution or a scalar-dominated solution. The critical points with 
k ^ arc cither curvature-scaling solutions or Milne universes. 

For critical points with k ^ 0, the curvature-scaling solution with k — —1 is the only 
stable solution. 

For critical points with A; = 0, the stabihty analysis shows that the fluid perturbations, 
curvature perturbations and scalar field perturbations all decouple from each other. The 
existence of a matter-scaling solution implies that there exists a scalar-dominated solution 
(when > 1/3). Only one of the two can be stable against fluid perturbations and 
that is the one with the fastest powerlaw. This implies that inflationary scalar-dominated 
solutions are stable against fluid perturbations since we assume that the barotropic fluid 
respects the strong energy condition. Only the inflationary scalar-dominated solutions are 
stable against curvature perturbations. 

Proper critical points are stable against scalar field perturbations in assisted infiation 
when all terms in the potential are positive. If some terms in the potential are negative there 
is no attractor. Since proper critical points have a higher powerlaw than the nonproper 
solutions, the critical point with the fastest expansion is the unique critical point that 
is stable against scalar field and fiuid perturbations. Furthermore, if the critical point is 
infiationary it is stable against curvature perturbations. Therefore in assisted inflation a 
fast expansion is dynamically favorable. 

For generalized assisted inflation the stability is more complicated. There are three 
main differences with assisted inflation. Firstly, nonproper critical points can be stable. 
Secondly, potentials with negative can have stable (nonproper) critical points, as long 
as the truncation of the negative terms can be done in a stable manner. Thirdly, there 
can exist unstable critical points with a higher powerlaw than the stable critical points. In 
that case the fastest expansion is not dynamically the most favorable. 

Despite the name generalized assisted inflation our models are more suited for present- 
day acceleration than early-universe inflation. Multiple exponential potentials either give 
rise to transient acceleration with too little e-foldings for early universe inflation or there 
exists an inflationary attractor without exit. 

Regarding present-day acceleration the shortcomings of our model are different. To 
avoid the cosmic coincidence problem, the period where the universe enters a stage of 
acceleration and where the energy density of the barotropic ffuid is comparable to that 
of the scalar fields has to be long enough. In our model such a period arises when the 
universe moves from an unstable matter-scaling solution towards the stable inflationary 
attractor [26]. During that transition there is a period with acceleration and comparable 
energy densities. But to make that period long enough one has to tune initial conditions. 
According to [27] this tuning becomes worse for multiple fields. A way to avoid this 
problem is by allowing interaction between the barotropic fluid and the scalar fields. Then 
there arises the possibility of inflationary attractors with nonzero energy density for the 
barotropic fluid [28]. 
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A Properties of the A-matrix 

Consider R hnearly independent vectors Their dual vectors that is, , = 5^ 
are uniquely determined when we demand that they all lie in Span{aa, a = 1 . . . R}. Then 
we have that 

[A-\, = {(5\P'). (96) 

The following R — \ vectors 

r = r- with a = l,...,i?-l, (97) 

lie in Spanjoa, a = 1 ... i? — 1} and are such that 

[A{R)-\, = 0\ 'P') . (98) 

Proof of the lemma 

The sum of all matrix elements of is given by 

ab ab<R a<R 

If we rewrite the /?'s in terms of the /3's via (fUTjl and use that , [3^) = 0, we find 

ab ab<R \H i H I ^^^^ 



a<R 

This can be rewritten as 

R 

IQR QR\ [E 

ab ab<R ; A- / 

which proves ()27|) since the second term on the RHS is positive. 



E[^"']'^^= E[^(i)"']'^''+7^7F^[E(/^^^'^)]'' (101) 



23 



Stability of truncations 

The condition for stability of truncating Yi is given by (j4T| . If we substitute 

Xi = -^JlY,aaAA{l)-\b, (102) 

^ ab 

in (pT|l the stabihty condition is given by 

Y,Aia{^aJb)>l. (103) 

a,b>l 

If we rewrite the /? in term of the /3 in ()l()4j) and then use the expressions ()9(ip98|l we get 

M 

a,b=l 

From which ()42|1 follows. 

Using ()42|) one proofs the stability conditions PHj) for truncating multiple Ya in an 
analogues way. 

Using the expressions ()96l98p one can proof with the same techniques as above that if 

M 

^[A-\b > for all b, (105) 

a=l 

then pSjl can never be satisfied. 

B Eigenvalues of didjW 

Here we will proof that didjW l^j,* has some negative values if some of the are negative. 

We use 0* to denote the value of the M — 1 scalar fields at the point where all first 
derivatives are zero, diW \(f,*= 0. For simplicity of notation we define the numbers r]a = 
Aq exp[— /{(eta, 0*)] and the M x M-matrix t] 

?7 = diag(?7i,?72, • • • , • (106) 

Similarly we define the vector fj as the vector with M components r/a. 

In this notation vanishing of the first derivatives implies the following M — 1 equations 

a,-f] = 0, (107) 

where we indicated the M — 1 columns of the matrix a as cJj. We further assume that the 
critical point we have in mind has p > 1/3 which implies that J2a Va > 0. 
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If the eigenvalues of the (M — 1) x (M — l)-matrix given by 



M 



didjW 



:i08) 



a=l 



are positive (p* is a stable minimum otherwise it is a saddle or a maximum. If all 77^ > it 
is straightforward to check that all eigenvalues are positive. 

Assume, after possible renumbering, that rji is negative, then define the M x M-matrix 
P as the matrix a with an extra column at the end with all entries of that column equal 
to 1 

/ ttll «12 • 
«21 ^22 



1 \ 

1 
1 

1/ 



(109) 



\ OLM\ 

The matrix a has maximal rank and therefore all columns are linearly independent. The 
same holds for the matrix /3. The reason is that the columns of a are all perpendicular to 
r/, whereas the last column with all components equal to 1 has the following innerproduct 
with rf 

((l,...,l),r^= J]r7,>0. (110) 



:iiii 



Hence (1, . . . , 1) must have a component in the rf-direction and it cannot be a linear com- 
bination of the columns of a. This means that (5 has maximal rank and since it is square 
we can define its inverse 

If we calculate the entries of M x M-matrix l3^r]i3 using the identity ()107|) we find that 
it takes the form 

We notice that the spectrum of (3'^'q(3 equals the spectrum of a^r]a plus the positive 
eigenvalue A = ^^^^a- If 13'^ has a negative eigenvalue so has a^r]a. Therefore it is 
enough to prove that there exists a vector %Ij such that ('?/'|/?^?7/3|'?/') < 0. It is easily 
checked that 

/ 1 \ 





'112 



V / 



is such a vector. 
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